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Abstract 

We explore several models which might be proposed to explain recent possible detections of 
high-energy (TeV) gamma rays in association with low-energy gamma-ray bursts (GRBs). Likely 
values (and/or upper limits) for the source energies in low- and high-energy gamma rays and 
hadrons are deduced for the burst sources associated with possible TeV gamma-ray detections by 
the Project GRAND array. Possible spectra for energetic gammas are deduced for three models: 1) 
inverse-Compton scattering of ambient photons from relativistic electrons; 2) proton-synchrotron 
emission; and 3) inelastic scattering of relativistic protons from ambient photons creating high- 
energy neutral pions, which decay into high-energy photons. These models rely on some basic 
assumptions about the GRB properties, e.g. that: the low- and high-energy gamma rays are 
produced at the same location; the time variability of the high-energy component can be estimated 
from the FWHM of the highest peak in the low-energy gamma ray light curve; and the variability- 
luminosity relation of Fenimore & Ramirez-Ruiz (2000) gives a reliable estimate of the redshifts 
of these bursts. We also explore the impact of each of these assumptions upon our models. We 
conclude that the energetic requirements are difficult to satisfy for any of these models unless, 
perhaps, either the photon beaming angle is much narrower for the high-energy component than for 
the low-energy GRB or the bursts occur at very low redshifts 0.01). Nevertheless, we find that 
the energetic requirements are most easily satisfied if TeV gamma rays are produced predominantly 
by inverse-Compton scattering with a magnetic field strength well below equipartition or by proton- 
synchrotron emission with a magnetic field strength near equipartition. 

PACS numbers: 98.70.Rz, 98.70.Sa, 95.55.Vj 
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I. INTRODUCTION 



Evidence has been accumulating for the arrival of ~ GeV-TeV gamma rays in coincidence 
with low-energy (~ MeV) gamma-ray bursts (GRBs). For example, EGRET detected seven 
GRBs which emitted high energy photons in the ~ 100 MeV to 18 GeV range . There 

have also been results suggestive of gamma rays beyond the GeV range J,0], although these 
were not claimed as firm detections. Evidence has also been reported for TeV emission in 
one burst out of 54 BATSE GRBs in the field of view of the Milagrito detector [6[ . In another 
paper Poirier et al. |7[ reported suggestive evidence for sub- TeV gamma rays arriving in 
coincidence with GRBs which occurred near zenith above the Gamma Ray Astrophysics 
at Notre Dame (GRAND) air shower array. In that experiment, most of the eight bursts 
analyzed were associated with at least some marginal excess (~ la) of muons including the 
event detected by Milagrito. One burst evidenced a possible detection at the 2.7a level. As 
shown in Poirier et al. if this detection is real, then the output in energetic gammas is 
likely to dominate the energetics of the burst. 

Although these data are not overwhelmingly convincing, they are at least suggestive that 
detectable energetic TeV gamma rays might be associated with low-energy gamma-ray bursts 
jjj. Moreover, these new detections, if real, can not be explained by a simple extrapolation 
of the BATSE spectrum ?|, particularly if intergalactic absorption is taken into account 
3- A new ~ TeV component in the GRB spectrum seems to be required. 

The present work is therefore an attempt to interpret these new possible detections in 
the context of three models, which might be proposed for the production of TeV gammas 
in a GRB. These are: 1) inverse-Compton scattering of ambient photons from relativistic 
electrons in the burst environment; 2) proton-synchrotron emission and 3) 



inelastic scattering of relativistic protons from ambient photons creating high-energy neutral 
pions, which decay into high-energy photons ^| . We here briefly outline the underlying 
physics and characteristic energetic gamma-ray spectra associated with each of these possible 
models. We then derive limits on the parameters of these models based upon the detection 
limits from the Project GRAND array. Based upon this, we conclude that it is difficult for 
any of these models to satisfy the energetic requirements unless the photon beaming angle is 
very narrow for the high-energy component. Of the models considered, the most likely are 
inverse-Compton scattering or proton-synchrotron emission. We note, however, that these 
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conclusions rely upon a few assumptions. For example, we have assumed that the low- and 
high-energy gamma rays are produced at the same location; that the time variability of 
the high-energy component can be estimated from the FWHM of the highest peak in the 
low-energy gamma ray light curve; and that the variability-luminosity relation of Fenimore 
& Ramirez-Ruiz (2000) gives a reliable estimate of the redshifts of these bursts. We explore 
the impact of each of these assumptions and find that, unless the bursts occur at very low 
redshifts 0.01), the energetic requirements remain difficult to satisfy. 



II. LOW-ENERGY GRB PROPERTIES 



The mystery of the astrophysical origin for low-energy gamma ray bursts (GRBs) has 
been with us for some time. As of yet there is no consensus explanation, although there 
is mounting evidence for an association with supernovae [16f. A likely scenario is a burst 

nnn 

environment involving collisions of an ultra relativistic e + — e plasma fireball jlTI . .18, .19]. 
These fireballs may produce low-energy gamma rays either by "internal" collisions of multiple 



shocks 



221. 



3, 



2l| , or by "external" collisions of a single shock with ambient interstellar material 



In either of these paradigms it is possible for very energetic gammas to be produced 
through inverse- Compton scattering of ambient photons off the relativistic electrons. Fur- 
thermore, it seems likely that baryons would be accelerated along with the pair plasma to 
very high energies 23 ] . Synchrotron emission from energeticprotons ]n ,12,13, 2^| 

or possibly hadro-production of pions in the burst environment 15j and subsequent 7r° 
gamma decay could also lead to an additional spectral component of very energetic gam- 
mas. In any case, it is at least plausible that energetic gammas could arrive in coincidence 
with a gamma-ray burst. This was the premise of Project GRANDs search for high-energy 
gammas arriving in coincidence with BATSE GRB observations. 

It is also possible, however, that low- and high-energy gamma-ray components are gen- 
erated in different regions or phases of a burst. This could lead to substantially different 
arrival times for each component. This was in fact the case for the 18 GeV photon observed 
by EGRET, which arrived ~ 4500 s after the low-energy emission had ended. 

Observations of TeV gammas could provide important clues as to the baryon loading, 
Lorentz factor, and ambient magnetic field of the relativistic fireball. Our goal in this paper 
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is to constrain the possible spectrum and source of energetic photons using the Project 
GRAND data. Hence, we restrict ourselves to considering high-energy gammas produced 
concurrently with low-energy gammas consistent with the search technique employed by 
Project GRAND. 



III. FITS TO OBSERVED GRB SPECTRA 

Table IJ summarizes some of the features of the BATSE and Project GRAND observations 
of the eight GRBs analyzed in Poirier et al. Q, where a detailed explanation of the Project 
GRAND results can be found. In the present paper, we will denote quantities in the frame 
of the observer by the superscript "06" . We will also use e 7 to denote the low energy GRB 
photons and distinguish them from the high-energy component, denoted E 1 . The observed 
BA TSE spectra are fit with a broken power law of the form 25[ 

, ob\ I ( ob,U^\ ( ,o6,MeV^-A -c c ob ^ c ob 



drf I / ob,MeV\-Ph -f c ob -> c ob 

I V C 7 / ' 7 — 7*' 

where ^ « 2, ^ « 1, and e° b b ~ 1 MeV is the break energy of the observed spectrum. 
Although these bursts are often better fit by using an exponential to join the two components, 
a broken power law is adequate for the present discussion. It maintains a simple analytic 
form for the equations, and as we shall see, the precise low energy form is almost irrelevant 
as long as a break energy exists. In what follows we will use values of a, Pi, and e ob b 
derived from optimum fits to the BATSE spectra for all events except the Milagrito event 
for which the BA TSE fluence was too weak to obtain a reliable spectral fit. The BATSE fit 
parameters corresponding to equation (0) are listed in Table UT1 The fit parameters for these 
bursts were provided at our request by M. S. Briggs at the Marshall Space Flight Center. 
We also include the variability time scale At for these bursts, which was estimated as the 
full width at half maximum (FWHM) of the brightest peak in each light curve. The light 
curves of GRBs typically show a wide range of timescale variability, so this choice may not 
be justified. In Will we explore the dependence of our conclusions on a wide range of values 
for At. 

In Table UTT) we give estimated redshifts for each burst. GRB 990123 is the only burst in 
this group for which an optical counterpart was detected. This burst, therefore, is known to 
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have occurred at a redshift of z = 1.6 126(1 . The redshifts for most of the remaining bursts 

n 

were estimated using the variability-luminosity relation of Fenimore & Ramirez- Ruiz [27] . 
This method relies on the apparent correlation between the time variability of a burst, which 
can be measured from its light curve, and the absolute luminosity of the burst, which we wish 
to infer. This provides us with a straightforward method for converting GRB observables 
into luminosities and redshifts |27j ]. 

Following Fenimore & Ramirez- Ruiz j^, we first fit a quadratic polynomial to the 
background in the non-burst portions of the BATSE 64 ms four channel data (i.e. DISSC 
data). Let hi be the binned background counts from this polynomial fit. If gi are the 
observed binned counts during the actual burst event, then the net count is q = — We 
then rebin the counts by dilating the time samples by Y — (1 + z)/(l + where z is the 



redshift we wish to estimate and Zb is a baseline redshift. Following [27J, we take Zb = 2. 
The new, dilated net counts, Cj, represent what the time history would look like at z = Zb- 
The variability is then defined to be the average mean-square of the variations in C\ relative 
to a smoothed time history, as 

T/ _ v -0.24 1 (C~ < C >) 2 — Bj 

where C p is the peak of the dilated net count during the burst and < C > is the count 
smoothed with a square-wave window with a length equal to 15% of the duration of the 
burst. The y -0 - 24 term corrects the variability for the energy- dependence of the time scale 
of a GRB [^J. The Bi term (dilated background counts in a sample) accounts for the 
Poisson noise. The sum is taken over the N samples that exceed the background by at least 
5cr. The estimated variabilities are listed in Table UTTl 

Based apon the flts of F — & Ramirez-Ruiz 0, we caa relate this v a ria bll , y V 
to the peak isotropic luminosity L 2 56 averaged over 256 ms in a specified energy range, e/ jP 
to e U)P (i.e. L 2 56 erg s -1 in the 50 to 300 keV band). This variability-luminosity relation is 

^25 6 /(4vr) = 1.9 x 10 51 y a86 erg s" 1 . (3) 

This peak luminosity depends upon the redshift, the observed spectral shape, and the ob- 
served peak photon flux P 2 56 (also averaged over 256 ms and over the same energy range) 
as 

L 256 = 4nD 2 z P 256 < e 7 > , (4) 



where D z is the co-moving distance and < e 7 > is the average photon energy in the luminosity 
bandpass per photon in the count bandpass. For this work we compute luminosities for an 
isotropic burst environment, such that Q = 4n, where Q is the unknown opening angle of 
the burst. If GRBs emit in a jet, our inferred luminosities and energies are diminished by 

The co-moving distance D z for a flat Qm + = 1 model is simply given by 



D z = jr [n M (i + z>) 3 + n A ]- 1/z dz> . (5) 

-"o Jo 

For the purposes of the present discussion, we will adopt the currently popular Q M = 0.3, 
Q A = 0.7, H = 72 km s" 1 Mpc" 1 model HHQ- 

From the observed photon spectrum 
the average photon energy in the luminosity bandpass per photon in the count bandpass is 

Hi 0[e 7 ] de 7 

where e\ and e u are the limits on the BA TSE energy range (~ 20 — 1500 keV). We can now 
iteratively solve equations (J2HSJ) until the estimate for z converges. 

This method of estimating the redshift of GRBs was found [jy] to be consistent with 
other estimates that rely upon an apparent relation between the luminosity and the time 



lag between hard and soft energy peaks, 
methods argues in favor of their reliability 



n he correlation between these two independent 
. However, in ^Vlll we explore the impact of 
systematically larger and smaller redshifts on our conclusions. 

Additionally, we can estimate the effective 4ir luminosity at the source in the BATSE 
energy band time- averaged over the full T90 interval, L 7 , by 

L 1=i ,Dljy^, d ^ + Z K (7) 

This is the luminosity estimate we use for the rest of the work presented in this paper. Both 
luminosity estimates are listed in Table ITTT1 for each burst. 

A. GRB 970417a 

GRB 970417a was the one burst (of 54 in the field of view) for which the Milagrito 
collaboration reported evidence of TeV emission during the duration of this burst within the 
BATSE error circle [fl. For this reason, we have included it in our analysis. Interestingly, 



this was a relatively weak BATSE burst with a fluence of 3.9 x 10~ 7 erg cm" 2 in all four 
BATSE energy channels (> 20 keV). Because this is such a weak low-energy burst, it is 
difficult to obtain reliable fits for the BATSE spectral parameters. For this reason, we 
have instead used the following average values of all bright bursts from Preece et al. j^: 
Pi = 1.0, (3h = 2.25, and = 225 keV. We can then use these average parameters and the 
observed flux to fix the normalization, a. The weak signal in the BATSE band also prevents 
us from using the variability-luminosity relation described above to determine the redshift 
of this burst. Instead we adopt z = 0.7 based upon the analysis of Totani [loT ]. 

IV. THE MODELS 

A. Inverse- Compt on Spectrum 

One possible source for energetic gamma rays is the inverse-Compton (IC) scattering of 
low-energy ambient photons by relativistic electrons. Indeed, such IC photons are thought 
to be the source of observed high-energy photons from active galactic nuclei such as Mk-421 



The inverse-Compton-scattering spectrum is generally written as 



dE^ / jq J dE e 



oc EZ 1 dE e —^ / dxf e . s {x) , (8) 



where d& e /dE e is the spectrum of electrons characteristic of a Fermi mechanism, 

rf$ A \E- a , if E< E b , 

A I (9) 



dE 



EbE- - 1 , if E > E h 



where A is a normalization constant and a k 2± 0.2 [35]. The flux of the synchrotron 
photons f e - s is 

{x-i"- 1 )/ 2 , below the cooling break, 
(10) 
x a l 2 , above the cooling break, 

where x cx E^jE 2 . Evaluating the integrals, we arrive at the simplified expression 

(d$ 7 / dEj)ic oc E' 5 , (11) 
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where the spectral index 5 is given by 

{a + l)/2 , for £ 7 < E cooi , 
5 = { (12) 

(a + 2)/2 , for £ 7 > E cool , 

where a ~ 2 is the electron spectral index and -E coo i is the energy scale at which the electron 
cooling time becomes comparable to the system lifetime. The inverse-Compton cooling time 
in the frame of the shock is 

where 7 e is the electron Lorentz factor and T = 300r3oo is the bulk Lorentz factor, both in 
the shock frame; L° 6 51 = L° b /10 51 erg s" 1 where L ob = L^/il + zf; and At = (l + z)r d /T 2 c is 
the variability timescale (in seconds) for the GRB, where r d is the radius at which low-energy 
gamma rays are emitted, both measured in the observer frame. 

The total fractional power radiated in inverse-Compton photons can be estimated from 
the ratio of the expansion time to the cooling time: 

f - = <T^k (14) 



_ 2 7ei#5l(l + *) 
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= 3.9 x lO- 2 1 " 7 ^V • (15) 

r 5 At 

1 300^^ 

The cooling frequency corresponds to fic — 1. The relation between the Compton-scattered 
photon energy and 7 e in the observer's frame is E° b ~ ^e° l ', where again e° b is the energy of 
the BATSE photons in the observer's frame. For simplicity we will take e ob = e° b b . Therefore, 
the cooling break energy can be written as 



plO ^2 £ o6,MeV 

E ^°- m (W + *)° GeV - (16) 

Relativistic Klein-Nishina (KN) corrections to the Compton spectrum may be important 
in the GRB environment over some range of energies. For completeness we include these by 
introducing a parameter T K n = E° b e° b {l + z) 2 /Y 2 m 2 e c A . When T K n 1) the KN effect is 
important. We adopt Tkn = 1 as the the lower boundary of the KN regime. This implies a 
lower limit to the observed gamma energy at which the KN effects should be considered, 

r 2 

#kn = 24 300 , M v GeV. (17) 



In the KN regime the emissivity of IC radiation per electron is independent of the electron 
energy and the cooling time becomes proportional to the electron energy [see e.g. Q. The 
emissivity is reduced by a factor of ~ i^kn compared with the classical IC formula (see Ecj. 
[T5]). and the observed photon energy is E° b ~ T^ e m e c 2 / '(1 + z). We can thus define the 
cooling frequency in the KN regime, by setting = /icT^n = 1- This gives 

j oh 

SS° 6 = 140 7 ' 5 . 1 M v GeV, (18) 

cool p 2 A7-^ o6 ' MeV ^2 V 7 

1 300^H e 7 6 / 

where we have used E° b ~ r7 e m e c 2 /(l + z) for the observed Compton-scattered photon 
energy. 

In the KN regime, the cooling is efficient in the lower photon energy range of E 1 < E™ x , 
i.e., contrary to the situation in the non-KN regime. In the KN regime, the photon index 
becomes 



for E 1 < E™ , 

a + 1, for E 1 > E™^ . 
To summarize, the IC spectrum can be modeled as 

^o^-( Q + 1 )/ 2 



(19) 



d$j\ A ic 
dEj lc ? 



{E° b ) > it E 7 < E coo \ , 

(-^cooi) 1//2 (E° b ) *■ + ^ , if -Ecooi < E 1 < E KN , 

7 ; (20) 

(-^cooi) 1 ^ 2 (E KN ) a ^ 2 1 (E° b ) , if -Ekn < E 1 < E™ x , 

^cooi) 1/2 (^kn)" 72 - 1 («) (EfY {a+l) , if E, > E™ , 
where all energies are assumed to be measured in GeV. 

For typical parameters, all three quantities, E coo \, Ekn, and E™ x are around 1-100 GeV. 
In any case, the spectrum beyond > 100 GeV must be steep with S ~ 3. The resultant 
IC gamma-ray spectrum should cut off above 

F 5/2 A ,l/2 

E° b {™ ax w Ef ' max = 4.2 x 10 5 1/A 300 GeV, (21) 



due to the cut-off in the spectrum of ultra-relativistic electrons at energies above E { 



ob,max 
e 



B. Spectrum of Proton-Synchrotron Gamma Rays 



It is generally believed that the expanding plasma contains at least some baryons. Indeed, 
some baryons within the jet are required to increase the burst duration and luminosity [cf. 

10 



371 ]. In the region where the electrons are accelerated, protons may also be accelerated 
up to ultra-high energies > 10 20 eV ^| producing a spectrum characteristic of a Fermi 
mechanism, (eq. [0]). 

The possibility of energetic protons producing ~Te V g ammas by synchrotron emission has 



been discussed in a number of papers 
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11 



12 



13 



2M . This mechanism has the desirable 



characteristic that the low-energy photons produced by electron-synchrotron emission and 
the high-ener gy p hotons from energetic protons can be produced simultaneously in the same 
environment j3Q]- 

In this model one assumes that there is a magnetic field present in the burst environment 
with approximate equipartition between the magnetic energy density B 2 /8tc and the total 
energy density, U. That is, 

g-UP, (22) 

where £g is a fraction of order unity. Following Totani [l3j we assume an optimally efficient 
proton-synchrotron environment in which U ~ U p ~ (m p /m e )U 7 , where U 1 is the photon 
energy density of BATSE gamma rays in the frame of the burst. This latter quantity can 
be related to the GRB luminosity utilizing L 1 = 47rr^T 2 cU 7 . We can then rewrite Eq. (|22|) 
in terms of the variables introduced in the previous section. This gives 

A/ 2 (Tpb \ 1/2(1 , \2 

B = lAx 10 ^ B [L ^ 51 J } L + Z) G . (23) 

r 3 At 

The photon energy from proton-synchrotron emission in the observer's frame is then 

, YYlehB 

jjjoo _ P 

7,p_s (1 + z)m p c 

= (E; b fehB(l + z) 
rripYc 5 

where Y p = E p /m p c 2 is the relativistic gamma factor of the protons in the frame of the 
fireball. This leads to the desired relation between the observed proton energy and the 
observed gamma energy, 

< = Cx(<£ v ) 1/2 , (25) 

where, 

F 2 At 1 / 2 

C = 1.6x 10 9 1M l30 ° At GeV. (26) 

£* 7 (^ b 5i) 1/4 (i + ^) 3/2 
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The final quantity needed to derive the energetic gamma spectrum is the cooling rate due 
to synchrotron emission. In the frame of the shock, the synchrotron cooling rate is 



1 _ 1 dE p 



^p— sync Ep dt 



Ae A B 2 E°\l + z) 
9m^c 7 T 



7.5 x 10- 2 ^ 4(L y 5 l)3 A y /2 + " )7/2 «^ v ) 1/2 b- 1 . (27) 



300 

The fractional energy loss to synchrotron photons then becomes 

TAt 



fp-sync(E p 



(l + z)t 



p— sync 

rpob 

ZP \f Rob < rnob 

pb,p — s 

i , if e* > e*^, 



(28) 



where 

= 7 ' 5 x 10 t W GeV ■ (29) 

Finally, we note that a proton flux <& p (E p ) will produce an energetic gamma flux of 

$ 7 = ^U{E P )% . (30) 

The observed high-energy photon spectrum can then be derived from equations Q, (|30j). 
and J251, 

= rf^ 
d£ 7 rfEp d£ 7 ' 1 J 



7. \ -, I f Z7 ob ,GeV\ 7 / 7T.o6,GeV\~( a + 1 )/ 2 ;f pob ^ rpob 

d^\ I i/W-J 1^7 ) ' Xt ^7 <jfc 7fe,P" 



to yield 

-1/2 

^7 j _ ^p-sync£<2-a y J \^7<>>P-s J V~ 7 / ' 7 ' — 7°.P-»' 

where 74P _sync is a normalization constant to be determined from observations. The break 
energy is 

Kp-s = 2-1 * 10- 3 3/2 ^°° At GeV . (33) 

The gamma-ray spectrum should cut off above E° b p ™ x p» 15r 30 o/(l + 2) TeV |l3j due to 
the cut-off in the spectrum of ultra-relativistic protons at energies above 

pB/2 a .1/2 

E ob,max = 2.1 x 10 11 — 7- ^0 GeV . (34) 
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C. Spectrum of Photo-Pion Gamma Rays 



Along with undergoing synchrotron radiation, protons accelerated in the burst environ- 
ment may collide with photons in the expanding fireball to produce secondary pions, which 
subsequently decay into high-energy gammas and neutrinos. This source of gamma rays 
seems unlikely due to the fact that it results from a secondary strong interaction and therefore 
has a small cross-section relative to electromagnetic interactions. Nevertheless an estimate 
of this spectrum is straightforward, so we include it here. Another alternative possibility 
might be pion production via proton-proton collisions [cf. |20j. However, the proton density 
in the frame of the shock must be small to ensure a low optical depth for gammas. Hence, 
p — 7 collisions are favored over p — p, although as we will show, even this preferred reaction 
places unreasonable energetic requirements on the burst environment. 

Following Waxman & Bahcall Q], the energy loss rate due to pion production is 

1 1 ( ™ d E«, m{ E)Er (35) 



tir 2r^ J Eo J e/2T P e2 de 

where Eq fh 0.15 GeV is the threshold for pion production. In the first integral, is the 
cross section for pion production due to a collision with a photon of energy e 7 in the rest 
frame of the proton, and £(E) is the average fractional energy lost to the pion. The second 
integral is over the low energy GRB spectrum, where </> 7 (e 7 ) is the photon flux in the frame 
of the proton. 

The evaluation of t n can be simplified |l^ by integrating the pion production cross section 
over the broken-power-law GRB spectrum (Eq. PQ) transformed back to the frame of the 
expanding plasma. Approximating the integral over the pion production cross section by the 
contribution from the peak of the A-resonance [as in 3| we deduce t n for a general spectral 



power law index flf 

1 dUty @peak£,peak A-Epeafc 

t-K £yb (1 ~t~ Pi) fcpc.ak 



x mm 



9F f: ^ ^ 



E 



peak 



s- 1 , (36) 



where <J pea k ~ 5 x 1CT 28 cm 2 is the A-resonance cross section, while E peak = 0.3 GeV and 
AE'peafc ~ 0.2 GeV are the energy and width of the resonance, respectively. The fractional 
energy lost at the peak is ^ pea k ~ 0.2, and U-y is the photon energy density of BATSE gamma 
rays in the frame of the fireball, as before. 
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As before we estimate the fractional power radiated as 
FAt 



fn 



(1 + Z% 

4.5x10-^ L^jl + zf j x " , if Ef < E* 



(37) 



7b ^oo^ 1 | xl , if Ef > £g >ff , 



where e° b b ~ 1 MeV is the break energy of the two power laws of the observed GRB spectrum. 
The last factor in equation (l3*7j) describes a break in the proton spectrum. In the observer 
frame this break energy is 15] 

E ° P l = 1-3 x 10 7 r h»,Mcv GeV ■ (38) 
[i + z) e jb 

Roughly half of the energy lost by the protons goes into 7r +, s, which quickly decay into 
neutrinos and positrons (through yU + s). In this work, we have ignored the effects of these 
decay products on the emerging gamma-ray spectrum. For neutrinos this is reasonable since 
there is very little chance of them interacting further. The positrons, however, may influence 
the gamma-ray spectrum through positron-synchrotron radiation |24| or pair annihilation. 

The other half of the energy lost by the protons goes into 7r°'s, which then decay into 
two photons. The mean pion energy is £, pe akE p . When the n° decays, the energy is shared 
equally among the photons. Hence, each gamma ray has an average energy 

E-y = ipeakEp/2 . (39) 

Now from equation (|3ip 



d%\ _ t ( 2 y- 



— \ AE, 



7 

1-Ph 



dE-y J ^ y £,peak J 

( -\ i ( Eiob,GevV ( T?ob,GcV\Ph-a-l . f ^ ^ Z7 



a;d^\ v - / * * > > — ( 4o) 

(1 + A)" 1 (E° b > G *v)- a , if e* > E$ >7r , 



where is a normalization constant and 



A, = 4.5 x 1(T 4 ^ — . (41) 

The observed break energy E°l n in the pion decay gamma spectrum is given by 



lb,: 

Kb, = %T^°t « 1-3 >< 10 6 -—%^v GeV . (42) 

(1 + z) 2 e lb ' 
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Above this energy the spectrum should obey the d$ 1 /dE 1 ~ E~ 2 of the protons and below 
this break energy, the exponent should be harder by one power, i.e. dQ^/dE^ ~ E~ x . As a 
practical matter, photons with energy as high as the break energy will not be observed, as 
they will be extinguished by pair production as described below. 

V. PHOTON, PROTON, AND ELECTRON LUMINOSITIES AT THE SOURCE 

From the above it is clear that the three models considered here imply different spectral 
shapes for the high-energy gamma component. Figure Q compares the initial spectra for 
all three mechanisms normalized to reproduce the Project GRAND observations for GRB 
971110 (as explained in Section 7). This would correspond to the unrealistic limit of no 
self or intergalactic absorption. Nevertheless, this illustrates the fact that these mechanisms 
have significantly different energetic requirements. Clearly, the most favorable energetically 
are the inverse-Compton and proton-synchrotron models. Figure El further illustrates this 
point by reproducing the source proton and electron spectra required by the various models, 
again normalized to reproduce the Project GRAND observations. Also included in this 
figure is the source electron spectrum required to only produce the observed BATSE data 
for this burst, following the electron-synchrotron model. The Project GRAND result, if it 
represents a real detection, requires a much higher flux of electrons to produce sufficient 
high-energy gamma rays through inverse-Compton scattering. This is probably a troubling 
requirement for the inverse-Compton model. 

To calculate these spectra we have used the normalization of the gamma-ray spectrum 
from the observed muon excess to determine the normalization of the associated ultra- 
relativistic proton and electron spectra. This is straightforward for the proton-synchrotron 
and photo-pion models, since the proton normalization appears explicitly in our final ex- 
pressions (equations [ID] and |H2j)- For the inverse-Compton model, we have followed Sari 
& Esin [38] to estimate the electron normalization from the following approximate relation 
between the low-energy electron-synchrotron spectrum, assumed to be observed by BATSE, 
and the inverse-Compton spectrum in the energy range of Project GRAND 



We take the electron number density n to be 1 cm 3 . The electron-synchrotron spectrum 




(43) 
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is assumed to have a form directly comparable to equation (J32J). We evaluate equation 



e°l to find A e ~ s ^ c . 



at £ 7 = E°l ol and e 7 

On the other hand, if one attributes |l4| the observed cosmic-ray excess above 10 20 eV 
to the energetic protons accelerated in GRBs, then an independent estimate of the proton 



normalization at the source can be obtained. Following Waxman [141 we note that the 
observed cosmic ray flux above 10 20 eV is ~ 3 x 10~ 21 cm -2 s _1 sr corresponding 



to an average universal number density of energetic cosmic rays of ucr ~ 10 30 cm 3 . If 

V20 



this density is due to GRBs then the number of protons with energies greater than 10 eV 
produced per GRB must be 



N(E P > 10 20 ) = Ucr/^tcr ~ 5 x 10 44 , 
where m 2 x 10 _10 (^/70) 3 Mpc -3 yr _1 is the cosmological GRB rate 



(44) 



40 



41 



421, and 



tcr ~ 3 x 10 8 yr is the lifetime of protons with E p > 10 20 eV. With our nominal E~ 2 spectral 
form, we deduce an absolute normalization of the proton spectrum emerging from an average 
GRB source of dN p /dE p = 5 x 10 55 Eq^ v GeV -1 . This is to be compared with the number 
of photons with E > 1 MeV emerging from a nominal GRB. If we assume that 10 53 erg 



7 ' 



is released in gammas above 1 MeV for an energetic photon spectrum of dQ~ f /dE 1 oc E, 
then the normalized spectrum for gammas above 1 MeV would be dN^/dE^ = 6 x lO 58 ^^ 
GeV -1 . Thus, one expects about 1200 gammas per proton from such a burst. 

Since both the proton-synchrotron and photo-pion models are based upon the same un- 
derlying baryon content in the relativistic plasma, it is instructive to summarize the relative 
efficiency of these two mechanisms for generating energetic gamma rays in GRBs. Let us 
consider a typical burst with f3h = 2 and a proton spectrum with a = 2. Then in the region 
E° b ~ 1 TeV, the fractional energy loss into these two mechanisms is 



A . 1 . 7X10 -^ 1+ ^ 4 



r 6 Ay- 

1 300 AAt 



TeV 



and 



p— sync 



mm 



6.7 x 10 2 — 7,51 



E. 



ob 



1/2 



F 4 At 1 / 2 

1 300 AAt 



TeV 



(45) 



(46) 



The factor of 1.7 x 10 7 suggests that proton-synchrotron emission will dominate over the 
photo-pion production for the parameters and gamma-ray energy considered. Only for 
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sufficiently small £b (~ 10 -4 ) does the efficiency of the proton-synchrotron mechanism begin 
to deviate from unity as ^ 4 . 

A direct comparison between inverse- Compton scattering and proton- synchrotron emis- 
sion is not possible since these two mechanisms rely upon different underlying energy sources, 
i.e. relativistic electrons for inverse-Compton scattering and relativistic protons for proton- 
synchrotron emission. Nevertheless, we note that both of these methods can be very efficient. 
In the region E° b ~ 1 TeV, fic ~ 1 and f p - sync ~ 1 for the range of parameters considered. 
This suggests that inverse-Compton scattering provides a very competitive mechanism for 
the generation of energetic gammas. Furthermore, it makes no requirement on the pre- 
existence of a magnetic field or baryon loading in the fireball plasma. 

Another important difference among all of these mechanisms is their associated cut-off 
energies. As noted above, there should be a cut-off for the proton-synchrotron spectrum 
around E° b JZT ~ 15 TeV, corresponding to a cut-off in the ultra-relativistic proton energies. 
The inverse-Compton spectrum has a much larger cut-off at around E° b ^ ax ~ 400 TeV, 
corresponding to a cut-off in the relativistic electron energies. For pion decay, however, the 
spectrum may extend all the way to 10 7 TeV ^j], but with a break at around 1300 TeV. 
These differences have a large effect on the implied total source luminosities of the bursts, 
as is apparent in Figure [T] 

VI. PAIR PRODUCTION OPTICAL DEPTH 

The spectra derived above must be corrected for two effects, both of which are due to pair 
production by energetic photons. First, within the burst environment, energetic gamma rays 
will interact with other photons to produce e + — e~ pairs. If this process is highly efficient, 
~ TeV gamma rays may not be able to escape from the burst. Even if some photons escape, 
this self-absorption will affect the implied source luminosity. 

The second effect is due to absorption along the line-of-sight from the burst environment. 
Here the energetic gamma rays interact with the intergalactic infrared and microwave back- 
grounds. This effect can cause a dramatic shift in the spectrum of energetic gammas in the 
TeV range depending upon the distance to the burst. 
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A. Internal Optical Depth from Pair Production 



A photon of energy E 1 interacts mainly with target photons of energy e 7 ~ 2m 2 e c i /E 1 
in the shock frame. We can approximate the cross section for pair-production as 3<7y/16, 
where a? = 6.6 x 10~ 25 cm 2 is the Thomson cross section. Then the optical depth can be 
approximated as 

3 D\ d^[e,/(l + z)] At 
r^nt ~ Y^T^e, — ^—^ , (47) 

where Atr/(1 + z) is the width of the emitting region as measured in the shock frame. 

This formula is similar to that of Waxman & Bahcall except that this form takes 

into account the spectral break in the low-energy GRB photons. This break is important 
since it implies that the optical depth is proportional to ~ for E° b < 2m 2 c A T 2 1 'e°j ~ 50 
GeV, but roughly constant for higher energies. The Waxman & Bahcall [l5j result is only 
valid in the lower energy range. The proper energy dependence is important because the 
internal optical depth can be of order unity. In the case of GRB 971110, the internal optical 
depth at 100 GeV is ~ 3 implying that some energetic gammas could emerge. The thin 
lines in Figure^show how the three source spectra normalized to fit GRB 971110 would be 
modified by internal absorption. The change in the energy dependence of the internal optical 
depth is apparent above about 50 GeV. For the remainder of the bursts, we were unable to 
place meaningful constraints on the internal optical depth due to large uncertainties in the 
fit parameters, particularly (3h and V . For these bursts, we have taken the very optimistic 
assumption of r 77iin f = 0. 

B. Intergalactic Optical Depth from Pair Production 

Another important constraint on the observed burst spectra comes from the absorption 
of photons via pair production during collisions with the intergalactic infrared background. 
In this work we use a calculated optical depth for intergalactic absorption based upon the 
standard formulation [e.g. Qj. We use a model for the luminosity evolution of background 



light from Totani, Yoshii, & Sato |43j. The dust emission component is calculated assuming 
a dust emission spectrum similar to that of the Solar neighborhood. The fraction of light 
absorbed by dust is adjusted to reproduce the observed far infrared background from COBE 
Q|. This method is summarized in Totani [ijj]. It is consistent with other optical-depth 
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calculations [cf . |9j, |45( . 

The resultant e + — e~ pairs may then further modify the original gamma-ray spectrum by 
providing a medium for some of the remaining gamma rays to undergo intergalactic inverse- 
Compton scattering. This process results in complicated showers of secondary electrons and 
gammas. These secondary gammas may be observed, but over a much longer time-scale, 
since much of this secondary light traverses a longer path length. The flux from these 
secondary gamma rays is probably below the detection threshold of current arrays such as 
Project GRAND. Nevertheless, this might be a noteworthy effect. 

Figure El shows final spectra in the observer's frame when both internal and intergalactic 
absorption are taken into account, assuming the burst is arriving from a redshift of z = 0.6 
(appropriate for GRB 971110). Even though the source spectra are vastly different, the 
observed high-energy gamma spectra are quite similar. Hence, the implied source energy 
requirement may be the only way to distinguish between the models. For illustration, Figure 
El also includes the observed BATSE spectrum for this burst, fit with both the Band (Eq. 
PQ) and electron-synchrotron models. 



VII. RESULTS 



In this work we use the muon observations of Project GRAND to fix the normalization 
(or upper limit) in each of the models described above for the various bursts analyzed. That 
is, the spectral shape is fixed from equations (jIUJ), (jH2J), or P0|). and the number of muons 
expected is then computed using 



r°° d$ (E ob ) 

N, = dAxT90x / dE°»P,(E; b ) £ j exp(-r 77 ) 



(48) 



where dA is the collecting area of the Project GRAND array (the effective area at the time 
of GRB 971110 was approximately 6.3 x 10 5 cm 2 ), E m i n is the primary gamma-ray detection 
threshold for Project GRAND (~ 10 GeV), and P M « 7.0 x 10 -5 (S^ eV ) 1 - 17 is the probability 
per primary for a muon to reach detection level at Project GRAND. This probability (valid 
for Ej > 3 GeV) was computed by Fasso & Poirier 46j using the Monte Carlo atmospheric 
absorption code, FLUKA. Here r 77 includes both the internal and intergalactic optical depth 
estimated for each burst as described above. For illustration, Table IIVI summarizes some 
estimates of the relative magnitudes of internal and intergalactic optical depths at two energy 



19 



scales, 100 GeV and 1 TeV. 

In practice, the integral in equation (}4*5j) is cut off at ~ 30 TeV since the optical depth is 
quite high for photons above this energy. The normalization constants for each model are 
then adjusted so that iV M agrees with the 2a upper limits set by Project GRAND, except for 
GRB 971110 where we used the observed mean value (see Table HJ). Here we adopt typical 
values for the degree of equipartition £b and the relativistic Lorentz factor T, namely £b ~ 1 
and T fa 300. Below we explore the dependence of our results on a broad range of these 
parameters. 

We can then use our normalized gamma spectra to estimate the total energy emitted in 
high-energy gammas at the source. For the inverse-Compton model, we can also estimate the 
energy emitted in electrons, noting again the characteristic spectrum of a Fermi mechanism 
given in equation Q. Similarly we can estimate the energy emitted in protons for the proton- 
synchrotron and photo-pion models. Implied energies for photon emission into An are given 
in Tables IVI[ and I VI II for the inverse-Compton, proton-synchrotron, and photo-pion 
models, respectively. We also list the required energies of the source protons or electrons, 
as appropriate. Our results assume a = 2 for the accelerated electron and proton spectra. 
For GRB 971110 we estimated the statistical uncertainties of our results using Monte Carlo 
techniques to explore the parameter space of each of the models assuming Gaussian error 
distributions. 

As evidenced by the large uncertainties, the models are not well constrained at present. 
Nevertheless, several points are worth noting from the tables. For the most significant 
possible detection (GRB 971110), the energetic requirements for the IC model (E^ ot,lc = 
(2.6 ± 8.6) x 10 55 erg) and the proton-synchrotron model (E^ ot ' p - sync = (3. ± 10.) x 10 55 
erg) are much less than that for a photo-pion mechanism (E^ ot,7T = (4. ± 27.) x 10 62 erg), as 
expected. 

The distinction between the inverse-Compton and proton-synchrotron mechanisms has an 
important consequence on the magnetic field of the GRB, specifically the degree of equipar- 
tition, £g. In the case of inverse-Compton scattering, the ratio of the inverse-Compton 
luminosity L 1 ^ to electron-synchrotron luminosity L sync should be equal to the ratio of the 
IC target photon energy density £7" 7 to the magnetic field energy density Ub since both 
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mechanisms arise from the same population of hot relativistic electrons. Thus, 




U B ' 



(49) 



L, 



'sync 



where Ub = £,bU and U is the total rest-frame energy density of the emission region. If we 
adopt the generally accepted view that the ~Me V gamma rays seen by BA TSE are caused 
by electron-synchrotron radiation then L sync = L T For GRB 971110, L 1 ^ jL 1 ~ 10 4 , which 
implies 



Since the energy density of the ~MeV BATSE gamma rays is likely to be higher than 
any other radiation source available for inverse-Compton scattering, we identify this energy 
density as U T This association has been implicit throughout this paper. On the other hand, 
since L 7 C ~ 10 4 L 7 the energy density of the ~TeV gamma rays must be greater than {7 7 by 
about 10 4 . We can then use 10 4 f/ 7 as a lower limit for the total rest-frame energy density U. 
Combining this with Eq. ()5fl)l we get the following upper limit on the degree of equipartition 



Note also that Eq. (}4T?|) holds only in the classical regime. In the KN regime £g is even smaller 
because inverse-Compton scattering becomes less efficient. Hence, the inverse-Compton 
model is at odds with models that propose GRBs as the source of ultra-high e nerg y cosmic 
rays, since those models require a magnetic field near equipartition, ~ 1 [of. Il4j . 

In Figures IH El and El we explore the dependence of our results upon a broad range 
of possible values for T, and At, respectively. These uncertainties were not formally 
included in our statistical estimates. Nevertheless, these figures support our qualitative 
conclusions. Specifically, the energetic requirements of all three models are in excess of that 
for the GRB for a broad range of parameters. Indeed most parameter variations in Figures 
0] -El only exacerbate this problem. However, the inverse-Compton and proton-synchrotron 
mechanisms are generally less sensitive to changes in T over the range 100 ^ T < 1000 or 
changes in At over the range 0.1 ^ At ^ 10. The sharp increase in the energy content of the 
protons and electrons at small V enters primarily through the internal optical depth. Figure 
El shows that the proton-synchrotron mechanism is highly efficient over a fairly broad range 



CbU < 10~ 4 f/. 



(50) 



6? < io 



8 



(51) 



of £ B (>0.01). 
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Regardless of the mechanism, if GRB 971110 is indeed a detection and our estimated 
redshift of z — 0.6 is valid, the implied energy in energetic (E 1 > 1 GeV) gamma rays is 
more than 100 times higher than the energy in the low-energy gamma rays. This renders 
the already challenging energetic requirements on the GRB source engine to be even more 
difficult. It may be possible to alleviate this difficulty if, perhaps, the photon beaming angle 
is much narrower for the high-energy component than for the low-energy GRB. Another 
possibility is that the redshift is much lower than our estimated value for this burst. In Figure 
[7|we explore the energetic requirements for GRB 971110 over a very broad range of redshifts 
from 0.005 to 3. This range is consistent with most currently measured GRB redshifts. This 
figure illustrates quite clearly the critical role that intergalactic absorption plays in driving 
up the energetic requirements and highlights the need for accurately determining the true 
redshifts of GRBs. 



VIII. CONCLUSION 

We have analyzed the eight GRBs discussed in Poirier et al. [t| , which occurred above the 
Project GRAND array. We have studied the implied energies in energetic (~TeV) gamma 
rays (and associated electrons and protons) of such GRBs in the context of three possi- 
ble mechanisms: inverse-Compton scattering, proton-synchrotron radiation, and photo-pion 
production. Our analysis suggest that all of these models face significant energetic require- 
ments. Gamma-ray production by either inverse-Compton scattering or proton-synchrotron 
radiation is probably the most efficient process. 

Although it can not be claimed that TeV gammas have unambiguously been detected in 
association with low-energy GRBs, we have argued that there is enough mounting evidence 
to warrant further study. Furthermore, we have shown that if TeV gammas continue to be 
observed, then they present some interesting dilemmas for GRB physics. In view of their 
potential as a probe of the GRB source environment, we argue that further efforts to measure 
energetic gamma rays in association with low-energy gamma-ray bursts are warranted. 
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10° 



— Inverse-Compton 




FIG. 1: Illustrative gamma-ray spectra for the three models discussed in the text. Spectra have 
been normalized to produce the muon excess observed by Project GRAND for GRB 971110. The 
curves begin at the detection threshold for Project GRAND and run to the respective cut-offs of 
each model. The thick lines are the raw source spectra. The thin lines illustrate the effects of 
internal pair-production optical depth on the source spectra. The change in energy dependence of 
the internal optical depth is apparent above about 50 GeV. 
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FIG. 2: Illustrative electron and proton spectra in the source frame for the three models discussed 
in the text. Spectra have been normalized to produce the muon excess observed by Project GRAND 
or to fit the observed BATSE spectrum for GRB 971110. 
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FIG. 3: Illustration of the effects of internal and intergalactic pair-production optical depth on the 
source spectra shown in Figure ^ This calculation assumes that the source (GRB 971110) is at 
a redshift of z = 0.6. The change in energy dependence of the internal optical depth is apparent 
above about 50 GeV. For illustration, this plot also includes the observed BATSE spectrum for 
this burst, fit with both the Band (Eq. pQ) and electron-synchrotron models. 
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FIG. 4: Plot of the energy escaping in gamma rays (accounting for r 77i j n t) and the source energy 
in protons or electrons for the inverse-Compton, proton-synchrotron, and photo-pion models as a 
function of the Lorentz boost of the GRB fireball. These curves correspond to spectra that were 
normalized to fit the observed muon excess for GRB 971110. The sharp increase in the energy 
content of the protons and electrons at small T enters primarily through the internal optical depth. 
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FIG. 5: Plot of the energy escaping in gamma rays (accounting for r 77i j n t) and the source energy in 
protons for the proton-synchrotron model as a function of the magnetic equipartition of the GRB 
fireball. These curves correspond to spectra that were normalized to fit the observed muon excess 
for GRB 971110. The offset between the photon and proton spectra for £b ~ 0.01 is due solely to 
absorption due to the internal optical depth, r 77i j nt . 
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FIG. 6: Plot of the energy escaping in gamma rays (accounting for r 77i j n t) and the source energy 
in protons or electrons for the inverse-Compton, proton-synchrotron, and photo-pion models as a 
function of the time variability of the GRB fireball. These curves correspond to spectra that were 
normalized to fit the observed muon excess for GRB 971110. 
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FIG. 7: Plot of the energy escaping in gamma rays (accounting for r 77i j n t) and the source energy 
in protons or electrons for the inverse-Compton, proton-synchrotron, and photo-pion models as 
a function of the redshift of the GRB fireball. These curves correspond to spectra that were 
normalized to fit the observed muon excess for GRB 971110. 
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TABLE I: Project GR AMD's Response to Selected BATSE Bursts 



GRB 


Trig 


T90 a 


RA a 


Dec a 


59 a 


N b 


971110 


6472 


195.2 


242 


50 


0.6 


467 ± 171 


990123 


7343 


62.5 


229 


42 


0.4 


< 75 


940526 


2994 


48.6 


132 


34 


1.7 


< 76 


980420 


6694 


39.9 


293 


27 


0.6 


< 133 


960428 


5450 


172.2 


304 


35 


1.0 


< 213 


980105 


6560 


36.8 


37 


52 


1.4 


< 107 


980301 


6619 


36.0 


148 


35 


1.3 


< 150 


970417a 


6188 


7.9 


290 c 


54 c 


0.5 C 


20 ± 17 



"Angles RA, Dec, and 59 in degrees and T90 in seconds. 
6 Upper limits are 2a confidence level. 

C RA, Dec, and 59 for GRB 970417a are based upon the Milagrito data. 



33 



TABLE II: Observed and Inferred GRB Properties in BA TSE Energy Range 





a 


ft 


Ph 




At 


GRB 


(cm- 2 s- 1 McV- 1 ) 






(MeV) 


(s) 


971110 


0.0751 ±0.0095 


1.02 ±0.04 


2.33 ± 0.11 


0.404 ± 0.041 


1.8 


990123 


1.93 ±0.14 


0.60 ±0.01 


3.11 ± 0.07 


1.29 ±0.07 


4.1 




< 0.217 


1.01 ±0.02 


> 3.2 


v 1 90 


1 2 


980420 


0.0815 ±0.0105 


0.18 ±0.08 


2.57 ±0.27 


0.553 ±0.033 


1.8 


960428 


0.0610 ±0.0086 


0.58 ±0.09 


2.49 ± 0.24 


0.433 ± 0.033 


0.7 


980105 


0.0332 ± 0.0081 


0.69 ±0.06 


2.82 ± 0.19 


0.286 ±0.031 


1.2 


980301 


0.0068 ± 0.0020 


0.54 ±0.20 


3.26 ± 1.41 


0.317 ±0.033 


2.3 


970417a 


0.0100 ±0.0061° 


1.0±0.5 i ' 


2.25±0.75 6 


0.250 ±0.150 i ' 


1.1 


"Spectral fits 


were not available for GRB 9704f 7a. Properties are based upon 


assumed redshift z ? 


a 0.7, 



observed BATSE flucnce, and average GRB spectral shape. 
b Average parameters for all bright GRBs considered in Preece et al. (2000). 
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TABLE III: Variability-Luminosity Redshift Estimates 





V 


td a 
X256 


z» 


Z/256 


C 


I 


GRB 




(photons cm -2 s -1 ) 




(ere: s" 


_1 ) 


(erg s _1 ) 


971110 


0.0204 


17.4 ±0.2 


0.6 ±0.2 


(6.6 ±5.2) 


x 10 51 


(3.8 ± 2.1) x 10 50 


990123 


0.0113 


16.6 ±0.2 


1.6 d 


(6.6 ± 1.0) 


x 10 52 


(7.8 ± 2.6) x 10 51 


940526 


0.0175 


14.4 ±0.3 


0.6 ±0.2 


(5.5 ±4.3) 


x 10 51 


(1.8 ± 0.7) x 10 50 


980420 


0.0267 


3.7 ±0.1 


1.1 ±0.3 


(7.2 ±4.8) 


x 10 51 


(1.0 ±0.4) x 10 51 


960428 


0.0628 


4.2 ±0.1 


1.6 ±0.5 


(1.9 ± 1.4) 


x 10 52 


(2.5 ± 1.1) x 10 51 


980105 


0.0671 


10.6 ±0.2 


1.2 ±0.4 


(2.4 ± 1.9) 


x 10 52 


(1.2 ±0.7) x 10 51 


980301 


0.0232 


0.9 ±0.1 


2.4 ±0.8 


(9.8 ±5.7) 


x 10 51 


(2.6 ± 5.0) x 10 51 


970417a 6 




0.7 ±0.1 


0.7 ±0. 2 f 


(4.3 ±3.1) 


x 10 50 


(0.2 ± 1.4) x 10 51 



Q Available at http://www.batse.msfc.nasa.gov/batse/ 

b The redshift uncertainties are crudely estimated as z/3 from the variability- luminosity relation. 

c We quote the isotropic luminosities. These must be multiplied by f2/(4-7r) to get the true luminosities. 

d Kulkarni et al. (1999) 

e GRB 970417a was too weak in the BATSE band to apply the variability-luminosity relation. 
■fTotani (2000). 
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TABLE IV: Internal & Intergalactic Optical Depths 



Internal Optical Depth 



Intergalactic Optical Depth 



GRB 



£ 7 = 100 GeV 



1 TeV 



E 1 = 100 GeV 



1 TeV 



971110 
990123 
940526 
980420 
960428 
980105 
980301 
970417a 



2.7 ±3.5 



2.9 ±3.6 



0.4 ±0.2 

2.8 
0.4 ±0.2 
1.3 ±0.8 
3.1 ±2.0 
1.6 ±1.2 
6.8 ±3.4 
0.5 ±0.3 



9.6 ±4.0 

31. 
9.6 ±4.0 

20. ±6. 
30. ± 10. 

22. ±8. 
44. ± 13. 

12. ±4. 



Targe uncertainties in the fit parameters made constraints on the internal optical depth uninformative in 
these cases. For these bursts, we have made the optimistic assumption r 77j i„i = 0. 
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TABLE V: Inferred Properties of Inverse-Compton Model 





Al c 

7 


COOL 


KN 


j-,KN,ob 
cool 


7pTot,IG a 

7 


T?Tot,lG a 

e 


GRB 


(cm- 2 s- 1 GeV- 1 ) 


(GeV) 


(GeV) 


(GeV) 


(ere) 


(ere;) 


971110 


4. ± 80. 


(1. ± 13.) x 10 7 


25. ± 7. 


75. ±48. 


(2.6 ±8.6) x 10 55 


(1. ± 15.) x 10 60 


990123 


< 6.2 x 10^2 


0.04 


2.8 


23. 


< 3.5 x 10 55 


< 2.6 x 10 57 


940526 


< 1.2 x 10~ 2 


14. 


7.3 


5.1 


< 2.5 x 10 54 


< 2.5 x 10 54 


980420 


< 4.2 x 10~ 2 


0.3 


9.8 


56. 


< 2.4 x 10 55 


< 1.6 x 10 56 


960428 


< 1.7 x 10~ 2 


0.004 


8.2 


380 


< 2.0 x 10 56 


< 2.4 x 10 59 


980105 


< 9.6 x 10~ 2 


0.04 


17. 


350 


< 3.9 x 10 55 


< 1.4 x 10 57 


980301 


< 2.1 x 10" 1 


0.05 


6.5 


73. 


< 3.4 x 10 56 


< 1.0 x 10 58 


970417a 


< 1.3 x 10~ 2 


8.6 


33. 


63. 


< 2.8 x 10 54 


< 2.8 x 10 54 



Here we estimate the total energy escaping in the gamma-ray component from each burst using the 
inverse-Compton model. We also estimate the total energy in the electron component at the source. 
"We quote the isotropic energies. These entries must be multiplied by Q,/{A-k) to get the true energies. 
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TABLE VI: Inferred Properties of Proton-Synchrotron Model 




Ap-sync 




rpTot,p—sync a 


rpTot,p—sync a 
P 


GRB 


(cuT 2 s _1 GeV -1 ) 


(GeV) 


(Vre"^ 




971110 


1. ± 22. 


(1. ± 22.) x 10 2 


(3 ± 10 ) x 10 55 


(1. ± 21.) x 10 58 


990123 


< 1.9 x 10~ 2 


6.2 x 10~ 5 


< 9 6 x 10 55 


< 9 6 x 10 55 


940526 


< 1.4 x 10~ 2 


1.3 x 10~ 2 


< 5 3 x 1 54 


< 5 3 x 1 54 


980420 


< 4.1 x 10~ 2 


8.9 x 10~ 4 


< 5.4 x 10 55 


< 5.4 x 10 55 


960428 


< 2.0 x 10~ 2 


5.8 x 10" 5 


< 2.7 x 10 56 


< 2.7 x 10 56 


980105 


< 3.8 x 10~ 2 


4.0 x 10~ 4 


< 5.6 x 10 55 


< 5.6 x 10 55 


980301 


< 9.3 x 10~ 2 


2.5 x 10~ 4 


< 7.0 x 10 56 


< 7.0 x 10 56 


970417a 


< 6.7 x 10~ 2 


3.0 x 10~ 2 


< 5.9 x 10 54 


< 5.9 x 10 54 



Here we estimate the total energy escaping in the gamma-ray component from each burst using the 
proton-synchrotron model. We also estimate the total energy in the proton component at the source. 
"We quote the isotropic energies. These entries must be multiplied by f2/(47r) to get the true energies. 
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TABLE VII: Inferred Properties of Photo-Pion Model 





Al 

p 




TpTot,TT a 
7 


TpTot.ir a 
P 


GRB 


(cm- 2 s- 1 GcV- 1 ) 


(10 5 GeV) 


(erg) 


(erg) 


971110 


(1. ± 13.) x 10 11 


13. ±4. 


(4. ± 27.) x 10 61 


(2. ± 27.) x 10 69 


990123 


< 1.3 x 10 20 


1.5 


< 1.2 x 10 74 


< 1.1 x 10 78 


940526 


< 1.8 x 10 13 


4.0 


< 2.0 x 10 65 


< 1.3 x 10 70 


980420 


< 7.8 x 10 15 


5.3 


< 2.5 x 10 69 


< 2.1 x 10 73 


960428 


< 3.9 x 10 18 


4.4 


< 5.3 x 10 73 


< 1.1 x 10 77 


980105 


< 2.3 x 10 17 


9.4 


< 1.8 x 10 71 


< 7.8 x 10 74 


980301 


< 3.1 x 10 27 


3.5 


< 3.4 x 10 81 


< 4.8 x 10 85 


970417a 


< 5.8 x 10 12 


18. 


< 3.2 x 10 64 


< 1.2 x 10 69 



Here we estimate the total energy escaping in the gamma-ray component from each burst using the 
photo-pion model. We also estimate the total energy in the proton component at the source. 
°We quote the isotropic energies. These entries must be multiplied by fi/(47r) to get the true energies. 
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